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ABSTRACT 

Inflation produces nearly Harrison-Zel'dovich scalar and tensor perturba- 
tion spectra which lead to anisotropy in the cosmic microwave background 
(CMB). The amplitudes and shapes of these spectra can be parametrized 
by Q s , r = Q\jQ 2 Sl rig and n T where Q s and Q\ are the scalar and ten- 
sor contributions to the square of the CMB quadrupole and rig and ny are 
the power-law spectral indices. Even if we restrict ourselves to information 
from angles greater than one third of a degree, three of these observables 
can be measured with some precision. The combination 130 1_ns Q| can be 
known to better than ±0.3%. The scalar index ns can be determined to 
better than ±0.02. The ratio r can be known to about ±0.1 for ns — 1 and 
slightly better for smaller ns- The precision with which ny can be measured 
depends weakly on ng and strongly on r. For ng — 1 ny can be determined 



with a precision of about ±0.056(1.5 + r)/r. A full-sky experiment with a 
20' beam using technology available today, similar to those being planned 
by several groups, can achieve the above precision. Good angular resolution 
is more important than high signal-to-noise ratio; for a given detector sen- 
sitivity and observing time a smaller beam provides more information than 
a larger beam. The uncertainties in ris and r are roughly proportional to 
the beam size. We briefly discuss the effects of uncertainty in the Hubble 
constant, baryon density, cosmological constant and ionization history. 



1 Introduction 



The detection of anisotropy of the Cosmic Microwave Background (CMB) 
by the COBE DMR [1] marks the beginning of a new era in observational 
cosmology. At least eight other experiments have subsequently made detec- 
tions at angular scales ranging from 0.5° to a few degrees [2]. The many 
possible sources of systematic error combined with the fact that most of the 
results are not much better than 3a detections means that care must be used 
in drawing conclusions from them [3]. However, it is encouraging that they 
are roughly consistent with each other and that several experiments have 
reproduced their results with repeated observations of the same area of the 
sky. The results are consistent with a nearly scale-invariant spectrum, pos- 
sibly with a feature at angular scales corresponding to the sound horizon at 
last-scattering (the "Doppler peak") [4, 5]. 

The current experiments, if not ruling out any theories with high confi- 
dence, are at the least providing strong constraints. For example, combined 
with large-scale structure data, the COBE DMR two-year data rule out the 
Primeval Isocurvature Baryon model [6] with 95% confidence or greater [7] 
and imply large bias factors for defect models [8]. 

In addition to testing theories, the microwave background also provides 
us with the opportunity to determine the parameters of a given theory. For 
example, in the cold dark matter (CDM) model the anisotropy depends on 
the amplitudes and slopes of the scalar and tensor spectra, the Hubble con- 
stant, Hq = 100/ikm/sec/Mpc, the baryon density, f^/i 2 , the cosmological 
constant in units of the critical density, Q\, and the redshift of reionization, 
zr. Bond et al. [9] have shown that, unfortunately, different choices of these 
parameters can lead to angular-power spectra which are indistinguishable by 
the current generation of experiments. This degeneracy makes it difficult 
to use CMB anisotropy experiments to determine cosmological parameters 
and hence the authors of [9] refer to the effect as "cosmic confusion". The 
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degeneracy has a positive effect as well. It limits the space of possible power- 
spectra, rendering the model testable despite its dependence on imprecisely 
known parameters [10, 11]. 

The long-term goals of the CMB observational community are much more 
aggressive than those of the current generation of experiments [12]. As a par- 
tial step toward those goals several groups in the United States and Europe 
are currently planning a "next generation" satellite experiment. Using detec- 
tor technology available today, such an experiment could sample the entire 
sky with a half-degree beam and in one year achieve a signal-to-noise ratio 
per beam-size pixel greater than one. Preliminary work by D. Spergel [13] 
indicates that such an experiment is sufficient for a significant lifting of the 
degeneracy of the angular-power spectrum pointed out by Bond et al. - at 
least in some regions of parameter space. 

Here we are interested in what a "next generation" satellite could tell us 
about inflation in particular. Besides resolving several cosmological puzzles, 
inflation is at the heart of the CDM scenario. In this picture, tensor, vector 
and scalar fluctuations in the metric are produced during an early epoch of 
rapid expansion driven by the vacuum energy of a scalar field [14]. The scalar 
perturbations grow via gravitational instability into the variety of structures 
we observe in the Universe today, as well as producing CMB anisotropy 
from about 10' scales up to the quadrupole. Vector perturbations decay 
with expansion and are of no phenomenological importance. The tensor 
perturbations today correspond to a stochastic background of gravity waves 
and also produce anisotropy in the microwave background at large angular 
scales (> 1°). 

Although there is no standard model of inflation, we can expect the spec- 
tra to have certain generic features. For inflation to occur, the dominant 
contribution to the energy density must be the vacuum energy of the scalar 
field. The kinetic energy is small in comparison, and hence the value of 
the scalar field changes slowly. Since the scalar field changes slowly while 
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the Universe is rapidly expanding, the perturbation spectra are nearly scale- 
invariant. Thus the power spectra are well-approximated by power laws with 
spectral indices close to the Harrison- Zel'dovich values. To be more precise, 
the primordial power spectra over the length scales of astrophysical interest 
at some time deep in the radiation-dominated era are well-approximated by 
the following power-laws in comoving wavenumber, k: 



with 0. As an example, the spectra from A</> 4 chaotic inflation 

[15] are best fit by n s = 0.94 and n T = -0.04. The fit is better than 0.5% 
in power from the quadrupole to the 10' scale. 

Since the two perturbation spectra are fit well by power laws, they can 
be characterized by four independent observables. We take them to be Q 2 S , 
r = Qt/Q"s, n s an d fiT-, where Q 2 S and Q\ are the expectation values of 
the scalar and tensor contributions to the square of the quadrupole. The 
quantities Q\ and Q\ should not be confused with the actual quadrupole 
moments on the sky. They are related in the same way that events from a 
random process are related to their parent distribution. If the perturbations 
are Gaussian (which is almost certainly the case for inflation), each of the 
five scalar (tensor) quadrupole moments on the sky is a single realization 
drawn from a Gaussian distribution with zero mean and variance Q\{Q%\ 

We wish to see how well these four inflationary observables can be deter- 
mined from a satellite experiment. The extent to which the other parameters 
(H , fl b h 2 , ...) can be determined as well remains an open question. Conceiv- 
ably, their confusing effects may be detrimental to the precision with which 
the inflationary observables can be determined. We must remember, though, 
that we have other sources of information on the cosmological parameters. 
Observations of light element abundances constrain VL^li 2 to be within the 
approximate range 0.009 to 0.022 [16] - a range which might well decrease 
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to ±10% in the next few years by the deuterium abundance measurements 
made in quasar absorption line systems [17]. The Hubble Space Telescope 
key project of calibrating Cepheids [18] and several physics-based methods 
[19] promise to make a definitive measurement of the Hubble constant in the 
near future to ±5%. Polarization of the CMB can provide constraints on 
ionization history [20]. Gravitational lens statistics constrain the value of 
a cosmological constant [21]. Thus in the following, we take specific values 
for the cosmological parameters and assume that they are perfectly known. 
Later we discuss how well they must be known in order that the uncertainty 
be negligible. 

It is worth pointing out that while the cosmological parameters may be 
determined by means other than CMB anisotropy, no observations are better- 
suited to determining the primordial spectra. Redshift surveys will continue 
to be plagued by theoretical uncertainties in the relationship of mass to light 
(the so-called bias), hampering determination of rig. Millisecond pulsars and 
space-based gravity wave detectors are probably not capable of detecting the 
very weak stochastic background of gravity waves expected from inflation - 
at least not in the near future [22, 23]. 

To simulate experiments, we need to assume particular values not only 
of cosmological parameters, but of Q 2 S , ns, tit and r as well. Below we focus 
on one case and then discuss how our results might change if the actual 
values are different. For the cosmological parameters we choose h = 0.5, 
Qa = 0, ilfe/i 2 = 0.0125 and the standard ionization history. For our theory 
of inflation we take the simplest model there is, chaotic inflation with a 
4> 4 potential. In addition to ns = 0.94 and rix = —0.04, chaotic inflation 
predicts r = 0.28. This is an example of a general rule for inflationary 
models called the consistency relation 1 r = —Itit- To choose Qs we note 

x This relationship, which is accurate to lowest order in (ns — 1) and nx, holds generally 
for slow-roll models; in some models (n — 1) ~ riy and the stronger relation, r = — 7(n— 1), 
also holds; see Ref. [24]. 
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that for ns = 1 and r = 0, the COBE DMR constrains the expectation value 
of the quadrupole to be Qs = 19.9 ±1.5 fiK [25]. Although the constraint 
would be slightly different for ns = 0.94 and r = 0.28, our only concern here 
is for rough agreement; we simply take Qs = 20 \xK. Chaotic inflation is an 
attractive choice for this study not only for its simplicity but also because r 
is twice as large as it has to be to ensure its detectability [26]. 

In section II we describe our calculation methods. We discuss the calcu- 
lation of the tensor and scalar angular-power spectra, and our modeling of 
experiments. In section III we show the results of attempts to recover Q s , r 
and ns from simulated experiments with varying beam sizes and signal-to- 
noise ratios. In section IV we see how well the consistency relation can be 
tested by attempting to recover all four of the observables. In section V we 
consider how our results would change if we had assumed different input val- 
ues of the inflationary observables and cosmological parameters. Particular 
attention is paid to the effect of a cosmological constant on the consistency 
relation. In section VI we briefly examine the effects of uncertainty in cos- 
mological parameters. Here we are interested in learning how well we have 
to know Hq, Q^Ii 2 , Q\, and reionization redshift, zr, so that our ignorance 
has a negligible effect on the determination of the inflationary observables. 

2 Calculation Methods 

The spherical harmonics provide a convenient basis for the expansion of 
CMB-temperature fluctuations: 

ST(6,<f>) = Y i aimY lm (e,<f>). (2) 

l,m 

Isotropy in the mean guarantees that (ai m a*, m ,) = Ci8u'5 mm i, where brackets 
indicate average over an ensemble of observers. It is the variance of the 
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multipoles that encodes information about the metric perturbations and C\ = 
(a 2 m ) is called the angular-power spectrum. (The expectation for the square 
of the quadrupole anisotropy is Q 2 = 5C 2 /47r.) Provided that the underlying 
perturbations are Gaussian, all predictions can be derived from the angular- 
power spectrum. 

For example, the expected value of the variance of temperature fluctua- 
tions from a given experiment is given by 



{ST2) = E *+± ClWl 



(3) 



where the window function W\ depends on the beam size and chopping strat- 
egy. For example, an experiment that measures the temperature difference 
between directions separated by angle 6 with beam size <r 6 has a window 
function W t = (1 - Pi(cos(9))e- 1 a i >. For a map made with a Gaussian beam, 
W, = e"' 2 < 




Figure 1: Tensor and scalar angular-power spectra for ns = 1 and rix = 0. 
Both spectra are arbitrarily normalized at C2 = 1. 
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We calculate the angular-power spectra by numerically evolving the pho- 
ton distribution function from deep in the radiation-dominated era until the 
present moment according to the first order general relativistic Boltzmann 
equation for radiative transfer. Details are given elsewhere [27]. The results 
of one calculation are shown in Fig. 1. 

To model the experiment, we assume that it creates a full-sky pixelized 
map of the CMB smoothed with a Gaussian beam with full-width at half- 
maximum #f w hm- The temperature of the i th pixel, 5Ti, has a contribution 
from the sky and from the instrument noise; STi = ST^ ky +5T* l01se . We further 
assume that the errors are uncorrelated and have uniform variance a 2 ix ; i.e., 
(ST™ ise 5T™ ise ) = a^Jij. For the moment we assume that 8T- ky = STp MB 
- an assumption we will soon discard. 

There are several different ways of describing the amount of noise in a 
map. The most straightforward way is to specify er p j x . Another way is to 
specify the signal-to-noise ratio per pixel, S/N. The observing time per pixel 
is inversely proportional to the pixel size, so both a p - lx and S/N depend on 
the pixel size. For definiteness, whenever referring to S/N or er p j x , we will 
take the pixel solid angle to be Q p - lx = #f w h m x #fwhm- To compare maps with 
different beam sizes, it is useful to have a measure of noise that is independent 
of fipix- For that purpose we use the weight per solid angle, w = (0-p ix fi p i x ) _1 . 

The error in each pixel, o- p j x , depends on the detector sensitivity s and 
the time spent observing each pixel, t p - lx ; a p i x = s/ \ft p \ x . The best detectors 
available today have sensitivities on the order of 200 //Ki/sec. With uniform 
full-sky coverage over the course of a year every 20' x 20' pixel could be 
observed for 85 seconds. Such a year of observing would result in a map 
with (7pj x = 22 \iK for fi p ; x = 20' x 20' - or a weight per solid angle of 
w = (7.5 fiK)~ 2 deg~ 2 . For comparison, the 2-year COBE maps have w ~ 
(400 /iK)~ 2 deg~ 2 . The considerable difference between these number is due 
to the ~ 70-fold improvement of detector sensitivities in the last 20 years. 

The signal used to calculate S/N is the rms of the temperature fluctu- 
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ations. The expected signal is given by Eq. 3. For the model we simulate 
here the expected signal with #f w h m = 20' is 92.5 fiK. Therefore the case of 
a pix = 22 jiK has S/N ~ 4. 

In our simulations, we never create a map. Instead we exploit the fact 
that the estimate of C[ which could be made from such a map, Cf Bt , would 
be xli+i distributed with mean (Cf st ) = C\ and variance 



(AQ) 2 = ({CT - C,)(CT - C V )) = ^ (Q + w-'e^f S w (4) 



(see appendix). 2 In the limit w = oo (cr p i x = 0), (AQ) 2 does not go to 
zero. This is because the finite sampling of events from a random process 
always leads to an uncertainty in the variance, called sampling variance, no 
matter how precisely each event is measured. The sampling variance for a 
Gaussian distribution is equal to twice the square of the variance divided by 
the number of samples. For each I there are 21 + 1 "samples" drawn from a 
Gaussian distribution of variance Q, hence the 2/(2/ + 1) factor in Eq. 4. In 
this limit of full-sky coverage, sampling variance is known as cosmic variance 



The signal at large / is reduced by the beam but the noise is not. If the 
beam profile is perfectly known, as is assumed here, one can take account of 
this diminution of signal by deconvolving the effect of the beam. The cost of 
doing so is the exponential factor in the noise term. 

Equation 4 can be rewritten in a more illuminating manner: 



This form is useful because l 2 Ci varies by less than an order of magnitude 

from I = 2 to / ~ 1000 for the models we consider. The cosmic variance 

2 Equation 4 does not include any error due to finite pixelization. At moderate to low 
S/N, these errors are unimportant if the pixel size is a few times smaller than the beam 
size. 



[28]. 
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term is proportional to \j\fl and dominates at small /. The noise term is 
proportional to I 3 / 2 at small / and for I > l/er 6 it increases exponentially. 

In Fig. 2 we show AC;/C; for four experiments with two different beam- 
sizes and two different values of w. One can see from the figure that, at 
constant w, the experiment with the better angular resolution is more precise 
at every value of /. The comparison at constant w is meaningful since these 
are experiments with the same detector sensitivity and observing time. 

From Eq. 4 it is easy to show that reducing the beam size at fixed 
detector sensitivity and observing time reduces AC; for every I, independent 
of Ci, as is evident in Fig. 2. At small values of I, AC; is near the cosmic 
variance limit and thus decreasing the beam decreases AC; only slightly. 
But for I > 1/(76 the reduction in AC; is dramatic. Thus there is much to 
be gained by reducing the size of the beam, even if this means reducing the 
signal-to-noise ratio to below unity. 

To this point we have assumed that 5T? MB = ST- ky . However, syn- 
chrotron and Bremstrahlung radiation, thermal emission from cold dust and 
unresolved extragalactic sources also contribute to the anisotropy of radiation 
at sub-millimeter to centimeter wavelengths at the angular scales of interest 
[29]. For this reason, a satellite experiment must make measurements over 
a range of wavelengths, so that the CMB component can be detected by its 
(hopefully) unique spectral dependence. Given 0- pix for a number of differ- 
ent frequencies, and guesses at the slopes and amplitudes of different con- 
taminating sources, one can estimate , the standard deviation in the 
determination of STp MB [30, 31]. Preliminary design studies by the MAP 
collaboration indicate that the foregrounds will degrade the noise level by 
a factor of two to three [32]. Therefore, we take into account the effect 
of foreground contamination by simply decreasing w from (7.5 /iK)~ 2 deg~ 2 
to (15 fiK)~ 2 deg~ 2 in one case and, to be conservative, (30 /iK)~ 2 deg~ 2 in 



9 



III 


, | ; , , 


V 1 1 1 1 1 h 1 1 1 1 ly 


— — - = 40' 






n on' 
° fwhm — flu 


1 1 
I 1 








cosmic variance 


III 




limit 

~~r— i— 


) 100 300 


300 


400 500 600 



1 

Figure 2: The precision of four different experiments. The dashed lines are 
for experiments with #f w h m = 40' and the dot-dashed lines are for experiments 
with f9f w hm = 20'. For each, the lower curve is for w = (15 /j,K)~ 2 deg~ 2 and 
the upper is for w = (30 nK)~ 2 deg~ 2 . 



another. 3 

Including the effects of foreground contamination makes the comparison 
of experiments with different beam sizes less straightforward. This is because 
the contamination may be more important at one angular scale than another 
and thus cause more degradation in a 40' experiment than in a 20' experi- 
ment. However, the smaller beam experiment will still always provide more 
information. As proof, we point out that it is always possible to synthesize a 
larger beam, after having done the experiment with a smaller beam. The in- 
teresting question becomes a quantitative one; how much better is a smaller 
beam? The answer depends on the spatial correlations of the contaminating 
sources, their frequency dependence and amplitudes and the frequency cov- 

3 The removal of foregrounds will not only increase the noise level but also introduce 
correlations in the noise from one pixel to the next. We ignore this effect [33]. 
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erage of the experiment. A complete investigation of this question is outside 
the scope of this paper. Here we only consider a very simple case of randomly 
distributed uniform point sources and an experiment that only measures at 
one frequency. For this case, the -\/N fluctuations in the number of point 
sources, N, in a given pixel imply (Tp™ b oc 1/Q pix . Thus for this case the 
comparison between different experiments should be made at equal weight 
per solid angle, just as for the no-foreground case. 

We simulate an experiment by drawing Cf st from the distribution in Eq. 
4 and then estimate Qs, ns and r from the "data" (the set of moments, 
{Cf st }) by finding the maximum of the likelihood function C(Qs, ns, r). The 
likelihood function is, up to a constant, independent of its arguments, de- 
fined as the probability density of measuring the set of moments Cf st given 
Ci(Q s ,n s ,r). 

C(Q s ,n s ,r) cx P(Cf st |Q(Q5, n s , r)) 



vf'- 2)/2 exp(-VV2) 

= 11; prr ttt^t — 1 where 

Q + w^e 12 ^ 2«'/ 2 r(n;/2) 

and rii = 21 + 1. The product on the right-hand side is simply that of the 
the xli+i distributions with mean C\ and variance given by Eq. 4 (see the 
appendix). 

To measure the certainty with which the observables can be determined 
we examine the distribution of the maxima from many simulations. An 
automated search, which uses the numerical technique of simulated annealing 
[34], finds the maximum for each set of simulated "data". Evaluation of the 
likelihood function on a fine grid of ns, Qs an d r shows that the maximum 
found by the automated procedure differs negligibly from the true maximum. 
Typically we perform 100 simulations which is sufficient to determine the 
standard deviation to better than 10%. 
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3 Results 
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Figure 3: The locations of the maxima from 100 simulations are shown pro- 
jected into the Qs — ns plane. The average values of ns and Qs, as well as 
their standard deviations, are indicated by the error bars. 

Figure 3 shows the maxima projected into the Qs — ns plane from 100 sim- 
ulations of an experiment with #f w h m = 20' and w = (15 /iK)~ 2 deg~ 2 . The 
average values of ns and Qs are equal to their input values to within 0.05% 
and hence there is no evidence for bias. The standard deviations in the values 
of ns and Qs are 0.016 and 0.74 /iK, respectively. We can conclude then 
that an experiment of this type can determine ns with Ans = 0.016 ± 0.001 
and Q s with AQ S = (0.74 ± 0.052) \iK. 

The strong correlation between Qs and ns, evident in Fig. 2, can be 
easily understood. It is due to the lever arm between I = 2 and those 
values of /, I*, for which C\ is measured most precisely. If C[* were the only 
moment measured then (Qs, ns) = (20 [iK, 0.94) would fit the data as well 
as (20 jiK (1* /2)(°- 94 ~ ns ^ 2 , ns) since each set of parameters results in nearly 
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Figure 4: The maxima of the likelihood functions for 100 simulations pro- 
jected into the Q* s (ns) — ns plane (see text). 



the same value of C;*. Keeping Q s oc (l*/2)~ ns ^ 2 as ns varies, causes the 
spectrum to "pivot" about I* and thus we call I* the pivot point of the data. 

The combination Q* s (ns) = o>Qs{l* /2) nS//2 is uncorrelated with ns- We 
choose the proportionality constant to be a = (/*/2)~ - 94 / 2 so that the mean 
value of Q* s equals the mean value of Qs- For the experiment under consid- 
eration, the pivot point is at I* = 210. Figure 4 shows the same maxima as 
in Fig. 3, but in the new coordinates ns,Q s (ns)- There is very little bias 
in the estimate of Q* s ; the average value of Q* s is 20.001 fiK. We find that 
AQ* s (n s ) = (0.0027 ±0.0002) \xK. 

For fixed C*, low Qs and high ns means that the fit has too little power 
at small /. To make up this deficit, r is overestimated. Thus there is also a 
correlation between r and ns (and hence Qs), a correlation which is evident 
in Fig. 5. We find Ar = 0.11 ± 0.008. There is very little bias; the average 
value of r is 0.29. 
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Figure 5: The maxima of the likelihood functions for 100 different simulations 
projected into the r — ns plane. 




Figure 6: The standard deviations of three of the four inflationary observables 
expected from experiments with three different beam sizes and two different 
weights per solid angle, w. Pentagons indicate Ar, squares Ans and triangles 
AQ* S /Q* S . The solid lines connect experiments with w = (15 fiK)~ 2 deg~ 2 . 
Along the dashed lines, w = (30 /j,K)~ 2 deg~ 2 . 
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We have analyzed the results from simulated experiments with a range 
of beam sizes and S/N ratios. The results are shown in Fig. 6 and Table 1. 
Each experiment has been simulated 100 times and thus the error bars are 
about the same size as the plotting symbols. 
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Table 1: The precision with which the three independent inflationary observ- 
ables can be measured. 

a) Units are (/iK)~ 2 deg~ 2 . 

b) analtyic 

c) numerical 



In order to understand our results analytically, it helps to take a greatly 
simplified view of the spectra. The approximate forms 

Z(Z + l)Cf = sciWY*- 1 ; 

l(l + l)C[ = §rCl{l/2)6{m-l) (7) 

capture the essence of the dependence of the exact spectra on Q s , ns, r 
and n T for 2 < / < 1000 and are sufficient for our purposes here. The 
detailed structure of peaks and troughs actually has little effect on how well 
the inflationary observables can be determined. 

The level of precision achieved by the simulated experiments can roughly 
be understood by considering measurements of the amplitude of the above 
spectra at three different angular scales centered around moments we shall 
call lc, Ib an d I a- The measurement at the smallest of the three angular 
scales is centered around lc and determines the amplitude of the scalar spec- 
trum. The multipole moment lc is defined to be that for which ACi/Ci is 
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a minimum. The measurement at the intermediate angular scale, centered 
around combined with the measurement at lc determines the slope. With 
the scalar spectrum thus completely determined, the value of r can be in- 
ferred from the amount of "excess power" in the large-angle measurement 
centered at I a- 

For simplicity, let us take the two smaller-scale window functions, Wf 
and W t B to be top hats with full widths lc and l B - The amplitude of the 
scalar spectrum can be known as well as we can know (ST C ). 



HQs) 2 „ A«ais» 
(Qs) 2 ~ 




(8) 

The 1 / sjlc factor is due to the reduction in variance caused by summing lc 
moments Q_; c / 2 to Ci + i c / 2 - Equation 8 accurately reproduces the numerical 
results for AQ* S /Q* S as can be seen in the Table. 

Measurement at a larger scale allows us to extract the slope. From Eq. 7 

n s -l 

->2\ 




m) - m) ( f ] (9) 

and therefore 

(Ans) 2 * (HI c /Ib))- z I ! ) + [ } ) (10) 




The analytic values of Ans m the Table are derived using Eq. 8 to find 
^ r2 y and then Eq. 10 with l B chose to minimize Ans- The analytic values 
consistently underestimate Ans- 

The experiment centered at I a can now be used to determine r. Knox 
and Turner [26] showed that a specific window function, W t A , centered at 
Ia — 55 is ideally suited for detecting the excess power due to gravity waves. 



16 



They used it to define an observable proportional to r: 

z - {STl) i fin 

\ J A /scalar 

where (£T|) sca i ar = Y,i ^^CfW^ is the anisotropy expected in experiment 
A from the scalar spectrum, which is fully determined by experiments B and 
C. If the proportionality constant ola is such that r = olaZa then 

(Ar) 2 = al(AZ A f 

I A{STj) \ 2 / A^rj) sc ,„ r \ 2 

= m 2 

Let us concentrate first on the case An,? = 0, the one considered in [26]. The 



((aA + r)AnsHl c /l A )) 2 . (12) 



last expression follows from the approximation of W A as a top hat centered 
at Ia with width Ia- Taking ola = 3 (since Cj 5 /C§ 5 ~ r/3, from Fig. 1) we 
find Ar = 0.055 (for r -C ola)- Knox and Turner found tmin = 0.14, where 
tmin is defined to be that value of r for which 95% of the time, r = can 
be ruled out with 95% confidence or greater. In order to compare results we 
must convert Ar to tmin- Since the tail of a Gaussian containing 5 per cent 
of the area is a distance 1.6tr from the maximum, tmin = 3.2 x Ar. Thus 
we find Tmin = 0.17, which compares fairly well with the exact result. 

The last expression for the second term in Eq. 12 follows from the fact 
that the uncertainty in the scalar contribution to Za is entirely due to Ans 
and the lever arm between I a — 55 and lc- The "analytic" results for Ar in 
the Table were calculated by substituting the numerical values of Ans into 
Eq. 12. 

4 Consistency Relation 

As mentioned above, inflation predicts a relationship between the tensor 
amplitude to scalar amplitude ratio and the shape of the tensor spectrum. 
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This relationship can be simply expressed in terms of the observables. To 
lowest order in the deviation from scale-invariance it is r = —Itit- If we 
could measure both r and rix with precision, this relationship would provide 
a powerful test of inflation [24]. Unfortunately, the only effect of the gravity 
waves that we can hope to detect in the near future is the increased anisotropy 
of the microwave background at large angular scales. This limited range of 
length scales over which the tensor spectrum has an observable influence 
makes the measurement of ny very difficult. 

For large r the tensor spectrum stands out more relative to the scalar 
spectrum, thereby decreasing Arix- For this reason, and because we are 
pessimistic about the prospects for determining n T well, we study the case 
r = 1, n T = —0.14. For the same reason, decreasing ns also decreases An T , 
but not as dramatically. However, the combination of large r and ns ^ 1 is 
strongly disfavored because it leaves insufficient power for structure formation 
on galactic scales. Therefore we choose to study ng = 1. 

In Fig. 7 the values of r*(%) and ny are shown which maximize the 
likelihood functions of 200 simulated experiments with #f w h m = 20', to = 
(15 yuK) _2 deg -2 . The combination r*{n T ) = r(P/2)" T+014 is analogous to 
Q*s( n s) m the scalar case. We find the pivot point to be at I* = 20 and 
An T = 0.14. 

The uncertainty in rix can be understood by imagining amplitude mea- 
surements at a pair of angular scales, exactly as was done for the scalar case. 
Rewriting Eq. 10 for ny 

(AM 2 - On(iVW)- 2 {(^If) + (^) 2 ) ■ (13) 

where Z D is defined analogously to Za- Equation 13 is correct only if we 
neglect the contributions to AZa and AZ D that are due to uncertainty in 
ns, since these contributions to the error in rix nearly cancel each other. The 
value of l D that minimizes An T is about Since Ia — 55, Id — 18. At 
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Figure 7: The maxima of the likelihood function of 200 simulations projected 
into the r*(n T ) — n T plane (see text). The "data" were synthesized for rig = 1, 
= 20 yuK, r = 1 and n T = —1/7. Unlike the previous cases n T is not 
constrained to obey the consistency equation r = — 7ny. 



this minimum, 

An T ~ (^^) 0-056. (14) 

From Fig. 1 we see that Cf 8 /Cf 8 ~ r/1.5. Therefore ckd ~ 1.5 and An^ = 
0.14 for r = 1, in agreement with the numerical result. 

Whether this precision allows for a test of inflation depends on the al- 
ternative hypotheses. For example, the hypothesis r = 1, tit = 1 is clearly 
distinguishable from the hypothesis r = 1, n T = —1/7. However n T = —r/7 
and n T = are not necessarily distinguishable for r < 1. We see that signif- 
icant deviations from the Harrison-Zel'dovich value, rix = 0, are necessary 
for the consistency relation to be falsified. 

Higher order corrections to the consistency relation are unlikely to be an 
important consideration in its testing. To second order in ns — 1 and tit, the 
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consistency relation is [35] 

n T = - T - (1 + 0.11r + 0.15(n s - 1)) . (15) 
Even for r = 2, the correction to the expected value of n T is only about 0.03. 

5 Dependence on Cosmological Parameters 

At fixed signal-to-noise ratio, the sensitivity of the above experiments to 
Qs and n,s is nearly independent of the actual values of h, flf,, Q\ and 715. 
However, ^a and rig do effect the sensitivities to r and n T through their effect 
on the shape of the scalar spectrum for I < 60. For example, decreasing ns 
at fixed r increases Cj jCf at I > 2 and hence improves sensitivity to r and 
riT- The sensitivity to ny also depends on r as shown in the previous section. 

The cosmological constant is the only cosmological parameter that signif- 
icantly affects the shape and amplitude of the scalar spectrum for / < 60. As 
the Universe expands and becomes cosmological constant dominated, the ex- 
pansion rate increases. The increased expansion rate causes the gravitational 
potential to decay which induces anisotropy through the Integrated Sachs- 
Wolfe (ISW) effect [36] , [37] . 4 The effect is largest for wavelengths that most 
recently entered the horizon and hence is largest at the quadrupole. For 
gravitational waves, the anisotropy in the radiation is all imprinted at the 
last-scattering surface [23] and hence a cosmological constant has little effect. 

The relationship between r and rix is due to a relationship between the 
primordial tensor and scalar spectra [24]. Therefore the dependence of Qs 
on 11a implies that the relationship between r and rix also depends on 
Oa- Defining fi A correction factors for the scalar and tensor quadrupoles, 

4 This has been called the late ISW effect by the authors of ref. [37] to distinguish it 
from the early ISW effect which occurred near the last-scattering surface as the Universe 
was in transition from radiation-domination to matter-domination. 
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fs(n A ) = (Qs(n A )/Qs(n A 
lows us to write the consistency relation for fl A ^ 0: 



0)) 2 and f T (Q A ) = (Q T (n A )/Q T (n A = 0)) 2 al- 



J T (n A ) 

-7 „ ,„ s tit- 



( 16 ) 

The correction factors as well as their ratio are shown in Fig. 8. They were 
calculated numerically using the Boltzmann codes described in references 
[27], extended to allow for a cosmological constant. The dependence of fr/fs 
on h and {l[,h 2 is much weaker. 
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Figure 8: Correction factors for the magnitude of the scalar and tensor 
quadrupoles, fs(Q A ) and f T (Q A ), as well as their ratio. 



6 Cosmic Confusion 

In addition to assuming particular values of the cosmological parameters, we 
also assumed zero uncertainty in their values. Here we briefly investigate 
how small these uncertainties must be so that they have a negligible effect 
on the precision with which we can determine the inflationary observables. 
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First we consider the confusing effect of uncertainty in Q b h 2 . Decreasing 
Q b h 2 decreases the peak at I ~ 200 (the "Doppler peak") and leaves the 
trough at / ~ 300 relatively unaffected. Thus a change in Q b h 2 can be 
compensated for by an increase in Qs, to fit the Doppler peak, and a decrease 
in n,s to fit the trough simultaneously. These shifts in Qs and ns create excess 
power at small I which is compensated for by an increase in r. 

To study these shifts in the inferred values of the inflationary observables 
we analyzed one set of simulated data several times, each time with a different 
assumption about the value of fl b . For #f w h m = 20', w = (15 fiK)~ 2 deg~ 2 we 
find the uncertainties in ns and r due to uncertainties in fl b h 2 to be 

(A%)o.». = °' 25 ^?. < 17 > 

Thus Q b h 2 must be known to better than 10% to be negligible for the mea- 
surement of r and better than 6% to be negligible for the determination of 
n s . 

To estimate the confusing effect of parameters other than Q b h 2 , we can 
use an empirical equation from Bond et al. [9] 

n w n s -0.281og(l + 0.8r)-0.515 ((1 - fl A ) h 2 ) 1/2 -0.000364 /2 + 0.26 (19) 

where Zr is the redshift of reionization (effectively for the standard ion- 
ization history). They claim that (for Q b h 2 = 0.0125) the angular-power 
spectrum is degenerate along surfaces of constant h, i.e. variations of pa- 
rameters at fixed h do not affect the shape of the angular-power spectrum 
significantly. To the extent that this is true, Eq. 19 implies 



{An s ) h = 0.52^/1 - Q A Ah = 0.23A/i, (20) 
(An s ) QA = ^j^- A^A = 0.47Afi A and (21) 

VI — l^A 
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{An s ) ZR = 0.13{z R /50)V 2 Az R /z R (22) 

where the left-most equalities hold for h = 0.8, Qa = 0.8. Therefore to get 
to Ans = 0.016 we must know the Hubble constant to better than 6% (14% 
for f^A = 0.8) and we must know £1a to better than ~ 7% (2% for 11a = 0.8, 
h = 0.8). 

A discussion of the prospects for precision measurements of the above 
parameters by means other than CMB anisotropy is beyond the scope of 
this paper. Clearly, given present uncertainties in cosmological parameters, 
the uncertainty in ns from any of the satellite experiments considered above 
would be "confusion-dominated". 

7 Discussion 

Our results should be compared to those of Hinshaw et al. [38]. They simu- 
lated full-sky maps of the CMB and then analyzed small patches by finding 
the maximum of the likelihood function for Qf, with all other parameters held 
fixed. They found, for S/N > 1, that sky coverage is more important than 
S/N and angular resolution. Angular resolution is their next most impor- 
tant factor; they found that Afl b decreases by a factor of 1.3 as the beam 
size decreases from 1° to 30'. Since they kept all other parameters fixed, the 
determination of flj, is effectively an amplitude determination and we can 
compare to our result for AQ* S /Q* S , for which, similarly, halving the beam 
size causes a factor of 1.3 decrease. 

This slight decrease should not, however, be used as an argument against 
the value of high resolution, the rewards of which are greater for the other 
observables. We find Ans and Ar are both proportional to #f w h m - Further- 
more, it should be emphasized, that even if it means reducing S/N to below 
unity, reducing the beam size (at fixed w) still results in dramatic improve- 
ment of our knowledge of the angular-power spectrum and the observables 
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studied here. 

Other considerations also argue for a small beam size. The extent to 
which the degeneracy pointed out by Bond et al. can be lifted probably 
depends critically on the beam size. Also, if Q + Qa < 1, the deviation of 
geodesies in an open Universe pushes the Doppler peak and all other features 
intrinsic to the last-scattering surface to smaller angles [39]. A lower limit 
on optimal beam size will probably come from constraints on the size of the 
telescope. 

We have seen that a full-sky map of the CMB at 20' resolution with 
S/N = 2 could achieve An s ~ 0.016, Ar ~ 0.11 and AQ*/Q* ~ 0.003. 
Such precision would allow for a critical test of inflation as the source of 
primordial perturbations. One particularly exciting prospect is the indirect 
detection of gravitational waves by determination of a non-zero r. While 
there is no generic inflationary prediction for the value of r, we note again 
that the simplest model gives r = 0.28 - a number significantly different from 
zero. Also, many models obey the relation r = —7(ns — 1) and a slight tilt 
(ns — 1 — 0.05) is helpful in fitting the large-scale structure data. Therefore 
r ~ 0.35 may be likely. Of course, there are also simple models that have 
negligibly small values of r. 

The effects of uncertainty in cosmological parameters (cosmic confusion) 
tempers our enthusiasm for the above precision. However, uncertainties in h 
and ilfe/i 2 are likely to decrease dramatically over the next few years. Also, 
precision measurement of multipole moments with I > 300 could possibly 
lift the degeneracy. Preliminary work by D. Spergel [13] suggests that this is 
indeed the case for a full-sky map with #f w h m = 0.5° and w = (10 fiK)~ 2 deg~ 2 . 
It would be interesting to see how the level of degeneracy changes as the beam 
size and weight per solid angle are varied. 
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Appendix: Calculation of ACi/Ci 

Let Xi be a Gaussian random variable with zero mean and variance a 2 . Then 
the sum 

n 

V^*?/* 2 (23) 

i 

is a random variable that has a \\ distribution. That is, it has probability 
density [40] 

y(n-2)/2 -V/2 

P(V)dV = .„ - . - dV. (24) 

v ' 2"/ 2 r(n/2) v ' 

To relate this result to the case of interest we make the identifications 

n = 21 + 1 and 

Xi = a^ AP (i = -> m = -I, i = n ^m = l) (25) 

where 

C AP =E^ m (^,^). (26) 
i=i 

There are two contributions to o^ lP - the signal convolved with the beam, 
a im e ~ l2a b/ 2 ^ and the noise, a°™ se . These two contributions are uncorrelated. 

i2 2 

The first, by definition of C/, has variance C\e a >> . The second has covariance 

2 

« se (« e )*) = 47T^^5 w (27) 
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which follows from applying the rules of error propagation to the definition 

MA 



of aj^ l AP and the assumption that the errors in the temperature of each pixel 



are uncorrelated with variance a 2 ix . Thus we make the further identification 

2 

a 2 = Qe- 12 ^ + 4tt-^P (28) 

since that is the total variance of a^ Ap . 

To estimate Q from C^ AP = £ m |a^ AP | 2 /(22 + 1) we must subtract off 
the expected noise contribution and correct for the finite width of the beam. 
Therefore 

C eBt = f C MAP _ 4n^-) . (29) 
Since the weight-per-solid angle is w = (o-p ix fipi x ) _1 = ^ 4 ^ gp ' x 

Cft = (pMAP _ w -lj e l*o* m (30) 



From this it follows that 



rieat _i_ ,,,-1 cpal 

V=W+l) C >+ W i e J (31) 



and therefore 



pander = p(v)^dcr 



n 



y(«-2)/2 e -v/2 



Cj + lo-V 2 '? 2-/2r( n /2) 



dCr (32) 



where n = 21 + 1. With the probability density for Cf st in hand, it is easy to 
show that (Cf st ) = Q and 

<(CT " C,)(CT - QO) = + ^- 1 ^ 2 ) 2 ( 33 ) 
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